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Abstract 



Tracking studies have indicated that for a lattice whose elements all have 
a single field multipole present, all having the same order k, the dynamic 
aperture approaches a non zero limit when k becomes very large. The dy- 
namic aperture and other properties of the lattice, as k becomes large, will 
be called the high multipole limit. It will be shown that the high multipole 
limit provides a reasonable estimate of the dynamic aperture of an acceler- 
ator, and the other properties of the high multipole limit found below are 
useful for understanding the stability of the accelerator. The high multipole 
limit is easily computed and it also provides an estimate of how much can 
be gained by correcting the lower field multipoles. The above results will be 
illustrated by tracking studies done with a simple one cell lattice, and with 
a RHIC lattice having six low beta insertions. 



Chapter 1 

Introduction 



Tracking studies have indicated that for a lattice whose elements all have 
a single field multipole present, all having the same order k, the dynamic 
aperture approaches a non-zero limit when k becomes very large. The dy- 
namic aperture and other properties of the lattice, as k becomes large, will 
be called the high multipole limit. It will be shown that the high multipole 
limit provides a reasonable estimate of the dynamic aperture of an acceler- 
ator, and the other properties of the high multipole limit found below are 
useful for understanding the stability of the accelerator. The high multipole 
limit is easily computed and it also provides an estimate of how much can 
be gained by correcting the lower field multipoles. 

The above results will be illustrated by tracking studies done with a sim- 
ple one cell lattice, and with a RHIC lattice having six low beta insertions. 
The properties of the high multipole limit that are demonstrated in these 
tracking studies can be used to answer the following kinds of questions about 
the dynamic aperture: 

1. Up to which order multipole does one have to correct to regain the 
aperture loss due to the field multipoles present. 

2. How much is the aperture loss due to the field multipoles present. 

The following are some properties of the high multipole limit discussed 
below which can be useful in understanding the stability of a given lattice: 

1. The stability boundary in the high multipole limit is dominated and 
determined by one set of elements in the lattice, and that is the set 
of elements which have the smallest value of -R//3 0,5 , where (3 is the 
linear beta function at the element and R is the multipole parameter 
that corresponds to the magnet radius. 
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2. Analytical results can be found for the stability boundary in the high 
multipole limit which can be used to estimate the loss in aperture due 
to the multipoles present in the lattice. 

3. The stability boundary in the high multipole limit does not depend on 
the strength of the multipoles or on the choice of the linear tunes, u x 
and v y . This leads to the suggestions that the stability boundary of 
a lattice is insensitive to the magnitude of the higher multipoles, and 
the linear tune needs to be chosen to avoid the resonances driven by 
the lower order multipoles. Higher and lower multipoles are defined 
below. 



Chapter 2 



The high multipole limit in 
2-dimensions 

This section will be devoted to establishing the basic rule for the high mul- 
tipole limit in 2-dimensional phase space. Consider a linear periodic lattice 
where each element of the lattice is perturbed by a single non-linear field 
multipole, and the order of this multipole, k, is the same in each element. 
The field multipole will produce a field in each element whose vertical com- 
ponent, By, in the median plane is given by 

B y = B b k x k , b k = b/R k (2.1) 

In computing the high multipole limit, we will be computing the dynamic 
aperture of this lattice for different values of k, and in particular for large 
values of k. R may vary from element to element, but does not vary with 
k. b may vary from element to element, and b may also vary with k but 
not by large factors. The dominant variation in b k with k is given by the 
1/R k factor. (More exactly, (Ab/b) 1 /^ approaches 1 for large enough k, 
where Ab is the largest change in b either from element to element or as 
a function of k.) For an actual accelerator, R may be chosen as the radius 
of the magnet coil, and the measured b = bkR k has to satisfy the above 
conditions in order to apply the high multipole limit results found below. 
Let xq, PxO be the initial particle coordinates. One may define the stability 
boundary to be a closed curve in xq, p x o such that for any choice of xo, p x o 
outside this boundary the particle motion for a given number of periods will 
be considered unstable. The definition of stability is discussed in section 6. 
The basic rule may now be stated as follows: 

Basic rule for the high multipole limit in 2 dimensional phase space. 
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For a particle moving through a linear periodic lattice in the 2 dimen- 
sional phase space of x, p x , in the presence of non linear field multipoles 
which have the form 

b k = h/R k (2.2) 

the stability boundary which encloses the stable area in xq, p x o, for a given 
number of periods and for large enough k, is given by 

t(X(),Pxo) = [R 2 /Px}min (2-3) 

e(x,p x ) = 7i i 2 + 2a x x p x + X pi 

m ; n is the minimum value of R 2 j ' (3 X in the elements of the lattice 
where bk is not zero, and j x , a x , Px are the linear parameters of the lattice. 
An analytical argument can be given which indicates what lies behind 



Eq. |2.3| . For very large k, the multipole field as given by Eq. 2.1 approaches 
zero when x in any element is smaller than the R value of that element. 
Thus for small enough x,e(x,p x ) becomes a constant of the motion. Some 
thought will then show that the element with the smallest value of R 2 / j3 x 
determines the largest emittance that is stable, which is given by Eq. 
A possible flaw in this argument is that for a given k, no matter how large, 
when x gets close enough to the stability boundary , the multipole field can 
become appreciably different from zero. 

The above basic rule will be justified below by doing a number of numer- 
ical tracking experiments, tracking particles through a number of different 
lattices. Most of the tracking experiments reported below are done with a 
simple one cell lattice. The results found will be further illustrated with 
results for a RHIC lattice with 6 low beta insertions. 



The simple one cell lattice 

The simple one cell lattice initially used in this study consists of a fo- 
cussing quadrupole, qf, and a defocussing quadrupole, q^, separated by drift 
spaces of equal length. The perturbing non-linear field multipole is initially 
placed in the middle of qf. The observation point for measuring the dy- 
namic aperture is initially chosen to be at the middle of the perturbing field 
multipole. This simple one cell lattice will be referred to as the simple one 
cell lattice. 

The perturbing field multipole is a point multipole which produces a 
vertical field on the median plane whose integrated strength, field times 
length, is given by 

B b x k /R k (2.4) 
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The parameters Bq, b, R are usually chosen to make the lattice resemble the 
RHIC lattice, with nonlinear effects of the same order as those seen in RHIC. 
To examine the high multipole limit one will be particularly interested in 
finding the dynamic aperture for large values of k. 

The transfer functions, that give the final particle coordinates for a given 
set of its initial coordinates for each element in the lattice, are given in 
section 7. For the reasons given there, the exact equations of motion are 
used in finding the transfer functions. The parameters of the quadrupole 
and drift spaces are initially chosen to produce the tune v x q = 0.1740. This 
tune was chosen to lie in a region free of all resonances up to the tenth order, 
and it lies between the 1/5 and 1/7 resonances. The lattice parameters are 
given in section 7. 

x s io vs k at qj for the simple lattice with a single bk at qf 

The basic rule for the hml (high multipole limit), and Eq. |2.3| will be 
illustrated by doing tracking runs using the simple one cell lattice. A single 
field multipole bk, as given by Eq. |2.4j , will be placed in the middle of the 
focusing quadrupole qf. The particle will be started with p X Q = 0, and x$ 
will be varied to find the largest xO that is stable for 100 periods, which will 
be denoted by x s iq. This will be done for different values of k, the order 



of the multipole. According to Eq. 2.3, since in this case p x o = 0, a x = 0, 



one should find for large enough k that x s iq = R. Here R was chosen as 



R = 0.04 m. The results are shown in Fig. 2.1 . 

One sees that at lower values of k, x s io varies rapidly and then levels out 
at a value close to x s iq = R = 0.04. The high multipole limit gives a good 
approximation for x s iq starting with the relatively low values of A; > 10 for 
non- linear multipoles of the order of those expected in RHIC. It will be seen 
below that this property helps to make the high multipole limit useful for 
estimating the dynamic aperture of an accelerator. 



The stability surface at qf for the simple lattice with a single bk at qf 

The stability boundary in the high multipole limit will now be found 
for the simple one cell lattice with just one bk at qf. To find the stability 
surface in xq, p X Q, one can search along different directions in xo,p X Q space 
to find the stability boundary in that direction. One may write the initial 
coordinates as 

xo = {Px € x0 ) ' 5 cos(air/2) 

Pxo = (e x o/p x )°- 5 sin(aiT/2) (2.5) 
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Figure 2.1: A plot of x s iq versus k for the simple one cell lattice with a 
single multipole of order k at qf. x s iq is the largest xq that is stable for 100 
periods. p X Q = p y o = yo = and 2 dimensional motion. In the figure xslO 
represents x s i . 



a gives the direction of search in xo, p x o and e x o is the initial linear emittance 
for this choice of Xq. p x o at qj where a x = 0, (3 X = 68.4497 m. To do 
the search e X Q will be increased until the motion becomes unstable for a 
particular choice of a. The value of e X Q that lies on the stability boundary 
for this search direction will be denoted by e xs iQ. According to the basic rule 
for the high multipole limit, Eq. 2J3, we should find that e xs iQ is constant 
for all directions at the value e xs io = R 2 /(3 X = 23.37482740. This would 
establish the basic rule for this lattice. This tracking study is done with a 
single field multipole, bf., at the middle of qj. In order to find the stability 
boundary for the high multipole limit, k is chosen at the large value of 
k = 1 • 10 6 . The results of this tracking study are shown in Table 2.1, where 
e xs io is shown clS cl function of a. 

One sees in Table 2.1 that e xs io for different directions is almost constant 



at the value 23.37482740 as predicted by the basic rule. However, there is a 
variation in e xs [Q of about .2%, which appears to show that the basic rule for 
high multipole limit, Eq. ^7J, is not exact but has a small error in it. This 
error is not important for the main results of this paper. This particular 
study shows that it is convenient to have a precise definition of stability such 
as is given in section 6. 

x s io at qd for the simple lattice with a single at qj 

Our next step will be to consider the stability limit at some other location 
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OL 




-1 


23.3893 


-0.75 


23.3749 


-0.5 


23.3756 


-0.25 


23.3748 


0. 


23.3747 


0.25 


23.386 


0.5 


23.4199 


0.75 


23.4024 


1. 


23.3893 



Table 2.1: The linear emittance, € xs iq, computed for different points, xq, PxO 
on the stability boundary in the hml, corresponding to different directions 
in xq, p X Q space, at the location of qf in the simple one cell lattice. The 
parameter a gives the different directions according to Eqs. p.5| . e xs io is in 
mm mrad. 

in the lattice, where (5 X is not at its maximum value, such as at qd where 
fi x is (3 X = 21.12345678. According to the basic rule for the high multipole 
limit for the simple one cell lattice, the stability limit in xo when p X Q = at 
Qfi Xsiof, is related to the stability limit at qd, x s id, at large enough k, by 

XslOd = (Pxd/Pxf) ' 5 XslOf (2.6) 

where (3 x f and (3 x d are the beta functions at qf and qd- This is illustrated 
in Table ^ where x s if and x s id are compared as a function of k, the order 
of the field multipole at qj. 



k 


Xslqd 


x slqf 


[Xslqd/ Xslqf] \J fixf / Pxd 


1000000 


0.02248406 


0.03999994 


1.0000187 


10000000 


0.0224841 


0.03999999 


1.0000193 


100000000 


0.02248411 


0.03999999 


1.0000198 



Table 2.2: A table showing how well x s i$ in the hml varies like f3 x ' 5 through 
the lattice. x s i a is the largest xo that is stable for 100 periods for the simple 
one cell lattice when p x0 = with a single multipole of order k at qf. x s iq at 
qf and qd are indicated by x s i f and x s iod, and j3 xq f and (3 xq d are the linear 
beta functions at qf and qd- All lengths are in meters. 
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One sees from Table 2.2 that the prediction of the basic rule for this case, 
as given by Eq. 2J3, is valid with an error of about 2 10~ 5 . The reason for 
varying k was to show that the error found was not due to k not being large 
enough as the basic rule for the high multipole limit holds only for large 
enough k values. The error in computing the beta functions also appears to 



be too small to account for the error found in Table 2.2 



x s io at qd and qf for the simple lattice with a single bk at qf and a bk at 
qd with a different R value 

We will now consider the case of the simple one cell lattice with two field 
multipoles present; one at qf and one at qd- The field multipole at qd will 



have an R- value (see Eq. 2A), R q d which is different from the R- value, R q f, 



of the multipole at qf. According to the basic rule for the hml, Eq. 2.2, 
when R q d = Rqf, the multipole at qf will dominate in determining the hml 
since [R 2 / P x ]mm will occur at qf where @ x has its maximum. In particular, 



Eq. gives x sU) w hen p s0 = as x s/0 / = Rqf, while x sl0 at q d is given by 
Eq. |2^ as x s /od = \/[Pxd/f3xf]Rqf- If one now reduces R q d, the multipole at 



qf will continue to dominate until R q d reaches the value y[(3 x d/P x f]Rqf, and 
then the multipole at qd will start to dominate, and x s iod will be given by 



Rqd while x s iof will be given by J[@ x f / Pxd\Rqd- These results are illustrated 



by the computed results given in Table 2.3 where x s ;o/ an d x s iod are shown 



as R q d/ Rqf is decreased from 1 to 0.2. For this lattice R q f is held constant 



at 0.04 m and ^[(3 xd /l3 x f] = 0.56209177. 

It will be suggested below that the high multipole limit provides a rea- 
sonable measure of the dynamic aperture. Assuming this to be so, then the 
basic rule for the hml states that in an accelerator the dynamic aperture is 
dominated by the magnet with the smallest R/(3®' 5 . For example, if an ac- 
celerator has 6 insertions with 6 crossing points, with similar magnets that 
are excited to give different (3 X at the crossing points, the magnet which 
has the largest (5 X will dominate and determine the dynamic aperture ||]. 
This magnet is located in the insertion with the smallest (5 X at the crossing 
point. The dynamic aperture at other locations in the lattice will scale like 
(/WArmax) ' 5 where P xmax is the largest X in the lattice. 

The dependence of x s io vs k at qf on the strength of field multipole, b 

In the previous tracking studies, the multipole present is assumed to 
have the form, bk = b/R k , where b is the integrated strength of the point 
multipole. b and R can vary from one element in the lattice to another. 
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^qd/ J^qf 


XslOqf 


7? j 

^qd 


XslOqd 


1. 


0.03999 


0.04 


0.02248 


0.9 


0.03999 


0.036 


0.02248 


0.8 


0.03999 


0.032 


0.02248 


0.7 


0.03999 


0.028 


0.02248 


0.6 


0.03999 


0.024 


0.02248 


0.55 


0.03914 


0.022 


0.02199 


0.5 


0.03558 


0.02 


0.01999 


0.4 


0.02846 


0.016 


0.01599 


0.3 


0.02134 


0.012 


0.01199 



Table 2.3: A table showing that the stability boundary in the hml is de- 
termined by the element in the lattice with the smallest value of R/(3®' 5 . 
Results shown are for the simple one cell lattice with a single multipole with 
the same order at qj and q^. R q j, R q( i and x s iof, x s ;od are the R and x s iq 
parameters at q* and qd- All lengths are in meters. 



The basic rule for the high multipole limit, Eq. 2J3, does not show any 
dependence on b. Tracking studies show that when b is increased, the high 
multipole limit is unchanged, but one must go to a larger value of k before 
the high multipole limit is reached. The studies also show that for large 
enough values of k, x s iq, the largest stable xo when p X Q = 0, goes 
This indicates that for large enough k, x s iq is insensitive to the size of b. 
For RHIC, the high multipole limit is reached for relatively low values of k, 
k ~ 10, and this result is insensitive to the size of the non-linear multipoles. 
Assuming the high multipole limit is a good measure of the dynamic aperture 
of an accelerator, then the dynamic aperture of an accelerator should be 
insensitive to the size of the higher multipoles, k larger than 10 for RHIC. 

Fig. 2^2 shows x s iq plotted against k for the simple one cell lattice with 
a single bt at qf. The two plots shown are for b = 0.024 m and for b = 0.24 
m. 

One sees that the same high multipole limit is reached at about k = 10 
for b = .024 and about k = 24 for b = 0.24. 

Fig. EOH shows thatx s /o ~ b^ 1 /*) by plotting x a i (Q.QM)/x lM (QM)/10 ( - 1 ^ 
against k. This ratio should approach 1 for large k if x s iq ~ b^^ k \ 

In the above it was shown that x s io ~ b( 1//fc ) which shows the dependence 
of x s iq on b for a fixed k. A more complete result which also shows the 
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Figure 2.2: A plot showing how x s ;o versus k depends on the strength, b, of 
the single multipole at qf for the simple one cell lattice, by comparing two 
cases where b differs by a factor of 10. In the figure xslO represents x s iq. 



dependence on k is 

x sl0 ~bW/k (2.7) 

The dependence of the high multipole limit stability boundary on the choice 
of linear tunes, v x and v y 

Tracking studies done with the simple one cell lattice indicate that the 
stability boundary in the high multipole limit does not depend on the choice 
of linear tunes, v x and v y . This leads to the suggestion, see section 8, that the 
linear tunes be chosen to avoid the resonances driven by the lower multipoles. 
The term lower multipoles is defined below. 
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Figure 2.3: A plot showing that x s iq depends on the multipole strength like 
x s io ~ b^ 1 /^ for large k, by comparing two cases where b differs by a factor 
of 10. In the figure, xslO represents x s iq. 



Chapter 3 

The high multipole limit in 4 
dimensions 



In 2 dimensional phase space motion, the stability boundary that encloses 
the stable area in xq, p x o in the high multipole limit is given by 

e(x ,Pxo) = [R 2 /Px]mm (3.1) 

One may ask what is the stability boundary for motion in 4 dimensional 
phase space in the high multipole limit. To answer this, one has to consider 
the motion of a particle moving in a lattice whose only nonlinear field in 
each element is that of a single multipole given by 

B y + %B X = -B b((x + iy)/R) k (3.2) 

where R may depend on s, and so may b although not by very large fac- 
tors. To find the high multipole limit stability boundary, one may use the 
argument given in section 2. For very large k, and for x 2 + y 2 not close to 
R 2 , Eq. |3^ shows that B y ~ B x ~ 0. Thus for x 2 + y 2 not close to R 2 , 

€x(x,p x ) = txO, e y(lJ,Py) = e y0 (3-3) 

where e x o, e y o are two constants and e x and e y are the linear emittance 
invariants. In addition for stable motion one has 

(x 2 + y 2 ) < R 2 (3.4) 



Equation ^4 can be restated as, using Eq. |3.3| , 

(Px e x0 + P y e y0 ) < R 2 (3.5) 
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Eq. 3.4 limits the range of e X Q, e y o for which the motion is stable. On the 
stability boundary , e x and e y are constant with the values of e x o and e y o 
respectively, and for each set of values of e x $ and e y o that are on the stability 
boundary, x 2 + y 2 must be less than or equal to R 2 for each element and at 
some point around the lattice x 2 + y 2 must be equal to R 2 . 

Eqs. 3J3 through |3.5| define the stability boundary in the high multipole 



limit. The stability boundary for the high multipole limit may be visualized 
as a curve in e^o, £y0 space as shown in Fig. |3.1| . For motion in 4 dimensional 
phase space, there does not appear to be a simple solution for the stability 
boundary in e x o, e y o space. This solution may depend on the particular 
form of f3 x (s), (3 y {s). There are, however, 3 points on the stability boundary 
for which one can find simple results. These are the 3 points for which 1. 
e x0 = 0, 2. e y0 = and 3. e x0 = e yQ . 



epyO vs epxO (mm. mrad.) 




epxO 



Figure 3.1: A plot showing the stability boundary in the hml for 4 dimen- 
sional motion, for the simple one cell lattice, as a plot of e y o versus e x o- 
epyO, epxO represent e y o, e X Q. 



For the e y o = point, Eqs. 13^ through give 

£x0 = [R 2 / Px]min, £y0 = 



(3.6) 



as the value of e x o on the stability boundary when e y o = 0. [R 2 /(3 x ] m in is 
the minimum value of R 2 /P x around the lattice. With this value of e X Q and 
e y = 0, one can show that (x 2 + y 2 ) is less than or equal to R 2 at every 
element in the lattice. Similarly for e x q = 0, one finds 



6,0 = [R 2 /Py\ 



j/Jmiii! 



CxO = 



(3.7) 
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When e x o = e v o, one finds 



xO — £j/0 — [R 2 / (Px + 0y)]mm (3-8) 



The three points on the stability boundary given by Eqs. |3.6| through 
provide a fairly good picture of the stability boundary. A fairly good approx- 
imation can be obtained by drawing two straight lines between the known 
three points. 

For accelerators which have insertion regions where X , (3 y have excep- 
tionally large values which occur in magnets which have the same R value, 
i?i ns , then Eqs. [Tl through |3.8] can written as 



£x0 — R 2 ns / flxma,x, £j/0 — 

Cj/0 = ^his/A/max, 6x0 = (3-9) 

t/)maK 



£x0 — CyO — Rfns/ (Px + A 



/^xmax is the largest j3 x in the lattice, and P yma . x and {fi x +{3 y ) ma x have similar 
meanings. 

Often, in tracking studies one does runs with initial values e x o = e y o and 
PxO = 0, p y o = 0, increasing xq until the motion becomes unstable for a given 
number of periods. The corresponding value of xq may be labeled x s iq. If 
one does tracking studies where all the bk present have the same k value, 
then for large enough k, the basic rule high multipole limit in 4 dimensions 



states that x s iq will approach a non-zero value given by Eq. 3.8 as 



XslO = ^J(Px0/R 2 )[R 2 /(f3x+Py)Un (3.10) 

where (3 x q and Rq are these parameters at the element where x s iq is mea- 



sured. For the insertion case described by Eq. 3.9, x s io is given by 



XslO 



{R /Rms)^/(3xo/{Px + A/W (3-H) 



In Fig. 3.1, e y o is plotted against e X Q showing the stability boundary in 
the initial emmitance space of e x o, e y Q. This curve was found by a tracking 
study using the simple one cell lattice with a single at qf. k was chosen 
at the large value of k = 1 10 6 , so that the curve is a good approximation 
of the high multipole limit. In this case the stability boundary in e^Oi £j/0 is 
almost a straight line. A simple approximation of the stability boundary is 
a straight line connecting the end points, the e x o = point and the e y o = 



point. According to Eqs. |3.6| through |3.8| this straight line is given by 
€xo/[R /Ac] min 

+ e yQ /[R 2 /f3 y ] min = l (3.12) 
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For the simple one cell lattice used here, where b^ is not zero only at the qf 
element, then [R 2 /f3 x ] min = R 2 / (3 xqf and [R 2 / (5y]mxn = R 2 /(3 yqf , R = 0.04 
m and (3xqf = 68.4497, (3yqf = 21.6265. One finds in this case that the 



three points on the stability boundary as given by Eqs. through agree 
with the tracking results with an error less than 5 10~ 3 . 
One may note that if [R 2 //3 x } min = [R 2 /P y ] 

m i n , which is true for RHIC 
because f3x = (3y at the low beta crossing points and tends to be true for 



proton colliders, then Eq. 3.12 gives 



CxO + e y0 = constant 



(3.13) 



on the stability boundary Tracking studies indicate that Eq. |3.13| is roughly 
true for RHIC. Assuming the high multipole limit provides a reasonable 
estimate of the dynamic aperture of the actual accelerator, then the above 
shows that the result that the total emittance, e x o + e y o, is roughly constant 
on the stability boundary is accidental in the sense that it depends on the 
properties of the beta functions at the crossing points. If these beta functions 
are not equal, then Eq. 3.12j will replace Eq. 3.13. 

Figure 3J2 illustrates the result given by Eq. 3.10 for x s iq for the case 
when e X Q = e y o, and p x o = p y o = in the high multipole limit. In fig. [3^ 
x s io is plotted against the multipole order, k. The results were found in a 
tracking study using the simple one cell lattice with a single at qf. For this 
lattice, with the parameters used, R = 0.04 m, and at qf, j3 x ,f3 y = 68.4497, 
21.6265 m. Eq. 3.10 then gives for x s i$ at very large k, x s iq = .0349 m. 
Fig. 3.2 shows that x s iq is approaching a value at large k near 0.0349m. At 
k = 20, x s ;o = .0348 m was found. 
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Figure 3.2: A plot of x s iq versus k for the simple one cell lattice with a single 
multipole of order k at qf, for 4 dimensional motion when p X Q = = p y o 
and e y o = exO. In the figure, xslO represents x s iq and is in meters. 




Chapter 4 



High multipole limit and the 
dynamic aperture 



The high multipole limit gives a result for the stability boundary that en- 
closes the stable area in xq, p x o- The goal of this section is to show that 
the high multipole limit gives a reasonable approximation for the stability 
boundary when all the field multipoles are present and the lower multipoles 
have been corrected. Conversely, the high multipole limit indicates how 
much may be gained by correcting the lower multipoles. The phrase lower 
multipoles will be more precisely defined below. Also one will see that when 
the nonlinear field multipoles are not too large, as is the case in RHIC, the 
high multipole limit provides a rough but useful estimate of the dynamic 
aperture. 

Motion in 2 dimensional phase space 

The above statements can be illustrated by the results of a tracking 
study in 2 dimensional phase space using the simple one cell lattice with 
non linear multipoles only at qt. To simulate the multipoles present in an 
accelerator, the point like nonlinear field at qj is given by 

By = B mathrmb{x/R) 2 {\ - (x/R) 50 )/{l - x/R) (4.1) 



Eq. ^4.1| gives a nonlinear field which contains all multipoles from k = 2 
to about k = 50, where all the multipoles decrease like 1/R k . With the 
parameters chosen as R = 0.04m, b = 0.024, this lattice resembles the 
RHIC accelerator without insertions. 



Fig. 4.1 shows the stability boundary in xq, p x o space as measured at qj. 



Two boundaries are shown; one is the high multipole limit, and the other 
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Figure 4.1: A plot of p X Q versus x Q showing the stability boundary for 2 
dimensional motion in the simple one cell lattice with all multipoles from 
k = 2 to k = 50 at qf. The stability boundary in the high multipole limit is 
also shown. In the figure, px0,x0 represent p x o, Po- 



boundary, that is contained inside the hml boundary, is the stability bound- 
ary when all the multipoles are present as given by Eq. |4.1| . According to 
the suggestions made at the beginning of this section, the difference between 
these two boundaries shows the loss in stable phase space due to the lower 
multipoles, and also how much phase space can be gained by correcting the 
lower multipoles. 

The term lower multipoles will be defined as follows. If one looks at 
Fig. 2A and Fig. |3.2| which plot x s io vs k for the two cases, e y o = and 
e y o = exO, one sees that x s iq gets close to the value given by the high 
multipole limit at about k = 10 for the simple one cell lattice. This value 
of k where x s iq gets close to the value given by the hml will be called fchmi- 
The lower multipoles are those multipoles for which k is less than fchmi- Let 
us now correct the lower multipoles from k = 2 to k = 9 giving the plot 



shown in Fig. 4.2 



Fig. 4.2 shows that by correcting 62 to 69, the stable phase space area has 
been increased so that it lies fairly close to the high multipole limit result, 
but still lies within the high multipole limit. If one corrects more multipoles 
past 69, the stability boundary will increase approaching the result for the 
high multipole limit. Results found using a RHIC lattice will be presented 
in section 6 which will also support the validity of the suggestions made in 
this section about the connection between the dynamic aperture and the 
high multipole limit. One may note that Fig. |4.1| shows that the loss in 



CHAPTER 4. HIGH MULTIPOLE LIMIT AND THE DYNAMIC APERTURE19 



P x0 pxO vs xO 
. 000£ i _ 



.0004 



hml 

N 



. 0002 



o.i04; 



02 



.0002 



. 



• * 

t 



02 f 



N 



.0004 



.000-6 



Figure 4.2: A plot of p X Q versus x Q showing the stability boundary for 2 
dimensional motion in the simple one cell lattice with multipoles from k = 10 
to k = 50 at qf. The lower multipoles, k = 2 to k = 9 have been corrected. 
The stability boundary in the high multipole limit is also shown. In the 
figure, px0,x0 represent p x o, xq. 



stable phase space due to the non linear multipoles used is about a factor 
of 2. About the same factor will be found for a RHIC lattice. Thus one can 
say that the hml provides a rough but useful estimate of dynamic aperture 
when the multipoles present are of the order of those expected in RHIC, 
overestimating the stable phase area in 2 dimensional phase space by about 
a factor of 2. 

Motion in 4 dimensional phase space 

The suggestions made at the beginning of this section can be illustrated 
by the results of a tracking study in 4 dimensional phase space using the 
simple one cell lattice with nonlinear multipoles only at qf. To simulate the 
multipoles present in an accelerator, the point like nonlinear field at qf is 
given by 

B y + iB x = B b((x + iy)/R) 2 (l - (z + iy)/R) 50 )/(l -(x + iy)/R) (4.2) 



Eq. ^4.2| gives a nonlinear field which contains all multipoles from k = 2 
to about k = 50, where all the multipoles decrease like 1/R k . With the 
parameters chosen as R = 0.04m, b = 0.024, this lattice resembles the 
RHIC accelerator without insertions. 

Fig. 4.3 shows the stability boundary in e x o, e y o space as measured at 
qf. Two surfaces are shown; one is the high multipole limit, and other 



CHAPTER 4. HIGH MULTIPOLE LIMIT AND THE DYNAMIC APERTURE20 



surface, that is contained within the high multipole limit surface, is the 



stability boundary when all the multipoles are present as given by Eq. 4,2 . 
According to the suggestions made at the beginning of this section, the 
difference between these two boundaries shows the loss in stable phase space 
due to the lower multipoles, and also how much phase space can be gained 
by correcting the lower multipoles. 
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Figure 4.3: A plot of e y o versus e X Q showing the stability boundary in e x o, 
e y o space for 4 dimensional motion in the simple one cell lattice when all the 
multipoles from k = 2 to k = 50 are present at qf. The stability boundary 
in the high multipole limit is also shown. In the figure, epy0,epx0 represent 

One may note that the high multipole limit stability boundary in e X Q, 
e y o space is a curve with zero thickness. However for the simple one cell 
lattice with all the multipoles present the e y o vs e x o stability curve has a 
certain smear or non-zero thickness. This is because for a choice of e x o, the 
e y o that lies on the stability boundary depends on the choice of xq, p x o, yo, 



p y Q. The curve shown in Fig. was obtained with p X Q = p y o = 0. 



Let us now correct the lower multipoles from k = 2 to k = 9 giving the 



plot shown in Fig. |4.4| . Fig. |4.4| shows that by correcting 62 to 69, the stable 
phase space area has been increased so that it lies fairly close to the high 
multipole limit result, but still lies within the high multipole limit boundary. 
If one corrects more multipoles past bg, the stability boundary will increase 
approaching the result for the high multipole limit. Results found using 
a RHIC lattice will be presented in section 5 which will also support the 
validity of the suggestions made in this section about the connection between 
the dynamic aperture and the high multipole limit. The results shown in 
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Figs. 4.3, and 4.4 also show how one can obtain misleading conclusions by 
looking at the results for just one direction in e x o, e yo space like the e x o = e y o 
direction. In this case the dynamic aperture actually became smaller for this 
direction when the k = 2 to 9 multipoles were corrected. 
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Figure 4.4: A plot of e y o versus e x o showing the stability boundary in e X Q, 
e y o space for 4 dimensional motion in the simple one cell lattice when the 
multipoles from k = 10 to k = 50 are present at qj. The multipoles from 
k = 2 to k = 9 have been corrected. The stability boundary in the high 
multipole limit is also shown. In the figure, epy0,epx0 represent e y o, e X Q. 



Chapter 5 

RHIC lattice results 



This section will illustrate the suggestions made in section 4 that the high 
multipole limit gives a reasonable estimate of the dynamic aperture when 
the lower multipoles are corrected by giving the results of tracking studies 
done with an early version of the RHIC lattice § . This lattice has random 
non linear multipoles in each element from k = 2 up to and including order 
k = 20. Skew multipoles are also present. The lattice has 6 insertions 
with 6 low beta crossing points at which (3 X = 6 m and /3 xmax = 236m, 
Pymax = 236m, {j3 x + /?j / ) m ax = 309 m. At qj in the normal cell, where x, 
Px, U, Py are observed, j3 x = 56 m, (3 y = 8.72m. 

The high multipole limit in RHIC 

In RHIC, the random multipoles in each element decrease roughly as 
R k . Because the multipoles are chosen randomly corresponding to given rms 
values ||, multipoles of different order or different k values, have different 
strengths or different b values. Elements with different R values are also 
present. The elements that have the smallest values of R 2 /f3 x or R 2 / (3 y , and 
which are the dominant elements, are in the insertions at the locations of 
flxmx and fiymx- The RHIC lattice differs from the simple one cell lattice 
in the presence of the chromaticity correcting sextupoles. In RHIC, besides 
the random multipoles whose rms values decrease like 1/R k , there is also 
a set of multipoles, the chromaticity correcting sextupoles, which do not 
fit into the 1/R k pattern of the random multipoles. Because of this, it is 
necessary to change the definition of the high multipole limit for RHIC. 
As described above, the high multipole limit is found by doing a series of 
tracking studies in which all the elements of the lattice have one multipole 
present of the same order, k, and the particle motion for very large k is 
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the high multipole limit. In the case of RHIC this procedure is changed 
in that in each tracking study the chromaticity correcting sextupoles are 
present as well as the random multipoles of order k. With this change in 
the definition of the high multipole limit, it is again suggested that the high 
multipole limit gives a reasonable estimate of the dynamic aperture when 
the lower multipoles are corrected. This will be illustrated by the following 
tracking results found using a RHIC lattice. This procedure for defining the 
high multipole limit can be used for any lattice for which there is another 
nonlinear field present as well as a set of multipoles that decrease like 1 / R k ; 
for example, one systematic multipole may be exceptionally large. 

This new definition of the high multipole limit changes the results found 
in sections 2 and 3 for the stability boundary in the high multipole limit. For 
very large k, and for x 2 + y 2 smaller than R 2 in each element, the particle 
motion is that of a particle in the presence of the chromaticity correcting 
sextupoles. In addition for stable motion one has 

{x 2 + y 2 )<R 2 (5.1) 

To find the stability boundary, one needs to know what is the maximum 
value of x 2 + y 2 for a given xq, p x o, yo, p y o. In section 3, this was given by 
the linear beta functions. In effect, one has to know what corresponds to 
the beta functions for the chromaticity correcting sextupoles for computing 
the maximum value of x 2 + y 2 . A semi-empirical solution of this problem 
will given below. One might notice one needs to answer this question only if 
one wants to have analytical results for points on the stability boundary in 
the high multipole limit like those given by Eqs. |3.6| through |3.8| . One can 
always find the stability boundary in the high multipole limit with tracking 
studies without much difficulty. 



Fig. |5.l| shows x s io plotted against k. In this study each element contains 
only one multipole of order k, and the multipoles in all the elements all have 
the same k value and the chromaticity correcting sextupoles are also present. 
Two curves are shown. For one curve e y o = 0, and for the second curve, 

x s iq is the largest xq that is stable for 500 turns when p y o = p X Q = 0. 



Fig. 5.1 is similar to Fig. 3.2 found for the simple one cell lattice and shows 
that x s io approaches a non zero limit as k becomes large. Using the results 
found by tracking runs one can find a result for computing x s io or y s io in the 
high multipole limit in RHIC. For the case when the chromaticity correcting 
sextupoles are absent, x s iq and y s iQ are given by Eq. |3^ as 

XslO 



(i?o/i?ins) \f Pxo/Pxma,x, £j/0 — 
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Figure 5.1: A plot of x s io versus k for RHIC. Each element has a single 
multipole of the same order, k. For one curve e y o = and p x o = p y o = 0. 
For the other curve e y o = e x o and p x o = p y o = 0. x s io is computed for 500 
turns. In the figure, xsl0,epy0,epx0 represent x s iq, e y o, e x Q. 



x sl0 — (Ro/ Rins)\J Pxol (Px + P\ 
(Ro/Rin 



VslO 



j/ymax; £y0 — ^xO 
Pyo/Pymux, ^x0 = 



(5.2) 



^j/max ; 

where P x o, (3 y o and Ro are these parameters at the element where x s iq is 
measured and for the insertion case described in section 3. To obtain a result 
that may be valid when chromaticity correcting sextupoles are present, we 



will replace Eqs. by 



XslO = fl (Ro I Rim ) \J PxO I Pxmax , CyO = 
%sl0 = H (Ro/ R\ns)\J fixo/ (Px + Py)max, £y0 
h (Ro/ Rins ) V PyO / P', 



£x0 



(5.3) 



VslO 



yvci&xi 



CxO = 



In Eq. p.3| , /i, /2, /3 can be found by using the results for xslO and ys£0 
found with tracking studies done for RHIC. This gives the results 



%sl0 — (Ro/Rms)y Pxo/Pxmax, £j/0 — 
X M = 0.707(Ro/Rins)^Pxo/(Px + Py) 
VslO = 



max; €y0 — ^x0 



(5.4) 



(Ro / Rms ) \J PyO I Pymax , £x0 — 



Eqs. 5.4 may be understood in the following way. The chromaticity 



correcting sextupoles in RHIC do not greatly distort the particle motion. 
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The stability surface in 2 dimensional phase space is a mildly distorted 
ellipse, as will be seen below. The main effect of the chromaticity correcting 
sextupoles is due to the coupling of the x and y motions, so that x 2 + y 2 
will grow from Xq + y^ by a factor which is found to be close to 1.414. Thus 



in Eqs. |TJ, the e y o = and the e x o = results are unchanged from those 
found when the chromaticity correcting sextupoles are absent, while in the 
£x0 = e^o case, the x, y coupling changes the result for xslO by the factor 
0.707. Although Eqs. |5.4| were found using tracking results for RHIC, they 
may be used for other proton storage rings when the chromaticity correcting 
sextupoles play about the same role in distorting the particle motion. Note 



that the factor of 0.707 in Eqs. [Lj is not based on analytical considerations, 
but was found though tracking studies with RHIC. 

Fig. 5^2 shows the stability boundary in the high multipole limit for 
motion in 2 dimensional phase space, e^o = 0. Tracking runs of 500 turns 
were used. The stability boundary is almost elliptical, showing that the 
chromaticity correcting sextupoles do not distort the motion very much. 

P x0 px0 vs xO 







9- 


..■oo-es 

.0004 


, 


****** 


Jiml 








0- 


. 0002 






N 


% 


4 














T 


\ 


03-0 


02-0 
0- 


01 

.0002 


. 


01 0. 


02 0. 


9 


\ 












y 








0- 


■ 0004 


— 




,* 





x0 



Figure 5.2: A plot showing the stability boundary for 500 turns in the high 
multipole limit for motion in 2 dimensional phase space for a RHIC lattice. 
In the figure, px0,xo represent p X Q, xq. 



Fig. |k3| shows the stability boundary in the high multipole limit for 
motion in 4 dimensional phase space, in e y o vs e x $ space. Tracking runs of 
500 turns were used with p X Q = p y Q = and with a single multipole with 
k = 20 present in each element. This surface has some thickness or smear 
which can be found by doing tracking runs with p X Q and p y o not 0. This 
curve is almost a straight line, and using Eqs. 5.3, it is described very well 



by Eq. 3.12 which in this case can be written as 

€xo/[R'ms/ Pxmx] ^yo/ [Rins / ftymx] = 1 



(5.5) 
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Figure 5.3: A plot showing the stability boundary in the high multipole 
limit for motion in 4 dimensional phase space for a RHIC lattice in e y o vs 
e x o space. Tracking runs of 500 turns were used with p X Q = p y o = and 
with a single multipole with k = 20 present in each element. In the figure, 
epy0,epx0 represent e^o, CxO which are in mm mrad. 



Motion in 2 dimensional phase space 

Tracking studies were done with a RHIC lattice to find the stability 
boundary for motion in 2 dimensional phase space, yo = PyO = 0. The 



results are shown in Fig. 5.4. Two curves are shown. The outer curve is 



the stability boundary in the high multipole limit as was shown in Fig. |5,2 . 
The inner curve is the path in x,p x for the last xo that was stable for 500 
turns as one increased xq with p X Q = and all the multipoles from k = 2 
to k = 20 are present. According to the suggestion being made here about 
the significance of the high multipole limit stability boundary, one would 
say that the lower multipoles have reduced the stable phase space by about 
36%. This loss in phase space can be recovered by correcting the lower 
multipoles, k less than about 10. Again, it is suggested here that the high 
multipole limit stability boundary indicates the stable phase space when the 
lower multipoles are corrected and it indicates the loss in phase space due 
to the lower multipoles. 



Fig. 5.5 is similar to Fig. 5.4 except that the multipoles from k = 2 
to k = 10 have been omitted. One sees that as the lower multipoles are 
corrected, the stability boundary approaches that of the high multipole limit. 
The loss in phase space has now been reduced to about 10%. Particle 
motions that came even closer to the hml and appeared to be stable for 
500 turns were seen and were rejected because of a rather large smear and 
scatter. 
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Figure 5.4: A plot showing the stability boundary for 500 turns for 2 di- 
mensional motion in RHIC with all the multipoles from k = 2 to k = 20 
present. The stability boundary in the high multipole limit is also shown. 
In the figure, px0,x0 represent p X Q, xq. 
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Figure 5.5: A plot showing the stability boundary for 500 turns for 2 dimen- 
sional motion in RHIC with the multipoles from k = 11 to k = 20 present. 
The multipoles for k = 2 to k = 10 have been corrected. The stability 
boundary in the high multipole limit is also shown. In the figure, px0,x0 
represent p x o, xq. 
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Motion in 4 dimensional phase space 

Fig. |5]6| shows the results of tracking studies done with a RHIC lattice 
to find the stability boundary for motion in 4 dimensional phase space. The 
stability boundary is shown by plotting e y o versus e^o for the case where 
PxO = PyO = 0. Two curves are shown. The outer boundary is the stability 
boundary in the high multipole limit. The inner boundary is the stability 
boundary when all the multipoles from k = 2 to k = 20 are present. The hml 
boundary was found by having only the k = 20 multipole present. Fig. [O^ 
shows a loss in 4 dimensional phase space of about 40% due to the presence 
of the lower multipoles. 
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Figure 5.6: A plot showing the stability boundary for 500 turns for 4 di- 
mensional motion for a RHIC lattice with all the multipoles from k = 2 to 
k = 20 present. The stability boundary in the high multipole limit is also 
shown. In the figure, epy0,epx0 represent e y o, e X Q are in mm-mrad. 



Fig. |5.7| shows the result when the multipoles from k = 2 to k = 10 are 
corrected. One sees that as one corrects the lower multipoles, the stability 
boundary approaches the high multipole limit stability boundary, and the 
loss in 4 dimensional phase space is reduced to about 10%. 
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Figure 5.7: A plot showing the stability boundary for 500 turns for 4 di- 
mensional motion for a RHIC lattice with the multipoles from k = 11 to 
k = 20 present. The multipoles for k = 2 to k = 10 have been corrected. 
The stability boundary in the high multipole limit is also shown. In the 
figure epy0,epx0 represent e y o, e x o and are in mm-mrad. 



Chapter 6 

Definition of stability 



In order to establish the properties of the high multipole limit, it is conve- 
nient to have a definition of stable motion which allows the stability bound- 
ary to be determined precisely In considering the motion of a particle in an 
accelerator, one might consider the particle motion for a certain number of 
periods to be stable if the particle motion stays within certain bounds, like 
those given by the vacuum tank, to be acceptable for the operation of the 
accelerator. Such a definition of stable motion, which allows a particular 
amount of acceptable growth in the particle motion, is not convenient for 
establishing the properties of the high multipole limit. A definition of sta- 
ble motion is given below which will precisely determine whether a particle 
motion for a given number of periods is stable. This definition may seem 
artificial. However, a good deal of numerical tracking shows that the stabil- 
ity limits found using this definition are usually close to the stability limits 
that would be acceptable for an accelerator. 

Consider the motion of a particle in a coordinate system which is based 
on a reference orbit where the independent coordinate is taken to be s, 
the distance along the reference orbit. The position of the particle is then 
described by x, y, and s, where x, y are the coordinates along two directions 
perpendicular to the reference orbit. The components of the momentum are 
then given by p x , p y and p s . If the energy of the particle is assumed to 
remain constant during the tracking, then p s can be computed from p s = 
{p 2 — pi — Py) ' 5 , p being the total momentum of the particle. The motion 
of the particle over a given number of periods will be said to be unstable 
if during the tracking of the particle over the given number of periods, 
p s becomes imaginary or p 2 + p 2 becomes larger than p 2 . p s becoming 
imaginary means that the formulation of the equations of motion based on 
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the given reference orbit has broken down because p s has changed sign, and 
the particle has reversed its direction along the reference orbit so that s is 
decreasing with time. 

The above definition of unstable motion over a given number of periods 
may seem artificial. However it has the advantage that the stability of a 
particular particle motion over a given number of periods can be precisely 
determined. Much experience with tracking also indicates that it is a use- 
ful definition and usually gives results that are close to the stability limits 
that would be acceptable for an accelerator. This definition of stability is 
convenient for establishing the above results for the high multipole limit. 
This precise definition of stability allows the stability limit to be calculated 
with great accuracy, and the tracking searches for the stability boundary 
can be automated. In order to use this definition of stability, one has to use 
the exact equations of motion. If one uses the approximations often used 
for large accelerators, where the radical (p 2 — p 2 — Py) ' 5 is expanded out 
assuming that p x /p and p y /p are much smaller than one, one will obtain 
invalid results as the expansion is not valid when p 2 — p 2 — p 2 is near zero 
and the radical is about to become imaginary. 

The definition of stability being proposed here has the following advan- 
tages: 

1. It avoids having to decide whether a particular particle motion is un- 
stable when some growth occurs and it is not obvious whether the 
growth is acceptable or not. 

2. With this definition of stability tracking searches for points on the 
stability boundary can be automated as it provides a simple test for 
stability. 

The results found in this paper do not depend on the choice of this definition 
of stability. The same results woud be found with any other reasonable 
definition of stability. 



Chapter 7 

Transfer functions for lattice 
elements 

In doing the tracking studies, one needs to know the transfer functions for 
each element of the lattice. The transfer functions allow one to compute the 
final coordinates of the particle from the initial coordinates for each element. 
As was indicated in section 6, the transfer functions have to satisfy the exact 
equations of motion in order to use the definition of stability given in section 
6. This can be accomplished by using the procedure |Q] of replacing a magnet 
with point magnets at the ends of the magnet separated by a drift space. 
By breaking the magnet up into pieces, one can approach the exact solution 
of the equations of motion by making the pieces smaller. One change in this 
procedure will be used here, which is that the reference orbit used will be 
made up of a series of smoothly joining straight lines and circular arcs ||. 
The transfer functions are then given in Ref. |5|. The circular arcs of the 
reference orbit are located at the dipoles in the lattice, and each arc has the 
curvature p which depends on the strength of the dipole. 1/p = at the 
quadrupoles and drift spaces. 

It is assumed that each magnet is broken up into a number of pieces. 
A magnet piece going from s = s\ to s = S2 and of length h = S2 — s\ is 
replaced by point magnets at the ends separated by a drift space of length 
h. In the following, q x = p x /p, q y = p y /p, q s = (1 - q 2 x - g^) ' 5 . 

A. Transfer functions for point magnets 

The transfer functions for a point magnet located at s = s\ is 

X2 = Xl, ?/2=2/l, 
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9x2 = fei + — -(l + (7.1) 

l /t sin 

9j/2 = g^i - -^-(l + xi/p)-^-5 ;E (xi,yi), 

/i is the length of the magnet piece, 9 = h/p. The field components B y and 
B x are assumed to depend only on x, y and do not change along the magnet, 
and that B s = 0. 

B. Transfer functions for drift spaces 

For a region along s in the lattice where l/p = for the reference orbit 

Qx2 = qxi, x 2 = xi + q x iLi2, 

q V 2 = q y i, V2 = yi + q y iL 12 , 

L\2 = (s 2 - si)/q sl (7.2) 

q S = (l-ql-ql)/ 1 ' 2 , 

L\2 is the path length between s\ and S2- 
For a region where l/p is not zero, 

q X 2 = qxi cos 9 + q s i sin 9, 
q S 2 = -q x isin9 + q sl cos9, 

9 = (s 2 - si)/p, 
x 2 = x 1 + {l + x 1 /p)2psm(9/2) (7.3) 
q x i cos 0/2 + q s i sin 0/2 
— g x i sin + q s i cos 
^12 = {l + x 1 /p)psm(9)/q s2 , 

q V 2 = q y i, 2/2 = 2/1 + <7j/i£i2 

C. Transfer functions for the simple one cell lattice 

This lattice has only point quadrupoles and drift spaces. For the transfer 
functions of the point quadrupoles one can use Eqs. 7A, replacing (h/2)B y 



and (h/2)B x by the integrated fields of the point magnet. For the drift 



spaces one can use Eqs. 1J2, The initial parameters that were used for the 
simple one cell lattice are the following: 

quadrupole integrated strength = 436.647 KG 

drift space length = 20m 

multipole field, b k = b / R k , b = 0.024, R = 0.04m 
Bp = 8400 KG. m 
B = 35 KG 



Chapter 8 



Longterm effects and the 
high multipole limit 

The stability boundary in the high multipole limit for 2 dimensional phase 
space does not appear to depend on nprd, the number of periods the particle 
is tracked. However many tracking studies have indicated that the stability 
boundary shrinks slowly the longer the particle is tracked. If one accepts 
the statement that the stability boundary in the high multipole limit is 
the boundary that is approached when the lower multipoles are corrected, 
then one can remove the apparent contradiction by the suggestion that the 
shrinking of the stability boundary, when nprd is increased, is due to the 
presence of the lower multipoles, and this effect can be reduced by correcting 
the lower multipoles. 

The following tracking study done with the simple one cell lattice sup- 
ports the previous statements. If one considers x s iq, the largest xO that is 
stable for a given number of periods when p X Q = 0, then one finds that x s iq 
decreases as nprd is increased. Using nprd = 10 2 and nprd = 10 , one finds 
the dxo/xo, the fractional decrease in x s io for these two values of nprd is 
dxo/xo = 0.033 when all the multipoles are present. If one corrects some of 
the lower multipoles by omitting the multipoles for k=2 to k=9, then one 
finds that (Ixq/xq is decreased by a factor of 6 to dxo/xQ = 0.005. 

Avoiding resonances of order 10 or higher 

It is sometimes suggested that in choosing the operating point for su- 
perconducting proton storage rings, one should avoid resonances of order 
10 or higher. A basis for this rule is provided by the high multipole limit. 
The range of the lower multipoles that reduce the dynamic aperture below 
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that given by the high multipole limit is given by the parameter k^mi de- 
fined in section 4. For RHIC, fc/, m ; is about khmi = 10. Since the important 
multipoles in affecting the dynamic aperture are the 10 lowest multipoles, 
it would seem desirable to avoid resonances up to 10 or higher which are 
the resonances driven by the 10 lowest multipoles in lowest order. If one 
would increase the strength of the non-linear multipoles by a factor of 10, 
thus raising khml to about khmi = 20, the above argument would suggest 
that one should avoid resonances up to order 20 or higher. 
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